eutad<-read.table("palmadulto.txt",header=TRUE, sep="\t")
head(eutad)
dist=matrix(NA, ncol=102, nrow=102)
for(i in 1:101)
{
  for(j in (i+1):102)
  {
    difx2=(eutad$gx[i]-eutad$gx[j])^2
    dify2=(eutad$gy[i]-eutad$gy[j])^2
    dist[i,j]<-sqrt(difx2 + dify2)
    dist[j,i]<-sqrt(difx2 + dify2)
  }
  
}
(nn<-apply(dist, 1, min, na.rm=TRUE))
(mnn<-mean(nn))

#exercicio
resultado<-rep(NA,1000) #criando 1000 NA
resultado[1]<-mnn #colocando mnn na primeira posicao

for(k in 2:1000) #repetindo 1000 vezes todas as operacoes
{
  xsim<-round(runif(102,0,320),1)
  ysim<-round(runif(102,0,320),1)
  dist.nova<-matrix(NA, ncol=102, nrow=102) 
  for(i in 1:101) # loop para diminuir a distancia entre o primeiro e segundo valor, segundo e terceiro, etc.
  {
    for(j in (i+1):102)
    {
    difx2=(xsim[i]-xsim[j])^2
    dify2=(ysim[i]-ysim[j])^2
    dist.nova[i,j]<-sqrt(difx2 + dify2)
    dist.nova[j,i]<-sqrt(difx2 + dify2)
  }
    }
    nn.nova<-apply(dist.nova, 1, min, na.rm=TRUE) #tirando os valores minimos
    mnn.nova<-mean(nn.nova)
    resultado[k]<-mnn.nova
}
hist(resultado)
abline(v=mnn,col="red")
valor.mnn<-sum(resultado<=mnn)
p<-valor.mnn/1000